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Abstract Recently the numerical simulations of the 
process of reionization of the universe at z > 6 have 
made a qualitative leap forward, reaching sufficient 
sizes and dynamic range to determine the character- 
istic scales of this process. This allowed making the 
first realistic predictions for a variety of observational 
signatures. We discuss recent results from large-scale 
radiative transfer and structure formation simulations 
on the observability of high-redshift Ly-a sources. We 
also briefly discuss the dependence of the characteristic 
scales and topology of the ionized and neutral patches 
on the reionization parameters. 
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Fig. 1. — Large-scale geometry of reionization and the 
large local variations in reionization history as seen at 
redshifted 21-cm line: (top) lg(5Ti)) and (bottom) 5Ti) 
smoothed with 3' compensated Gaussian beam. 




Fig. 2. — Ionization maps (blue neutral, red ionized) of 
selected 100 h~^Mpc box simulations for high (bottom 
panels) and low source efficiencies (top) and with (right 
panels) and without (left panels) sub-grid clumping. 
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1 Introduction 



Fig. 3. — The reionization history of a high density 
peak. The images are centered on the most massive (at 
z = 6) halo in our computational volume and are of 
size 100/i~^Mpc to the side. The snapshots are (left 
to right): z = 12.9, z = 9.0, z = 8.0, and z = 6.6. 
The underlying cosmological density field (dark, green) 
is superimposed with the ionized fraction (light, orange) 
and the ionizing sources (dark, blue dots). 
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Fig. 4. — Ly-a sources at redshift z = 6.0: Transmis- 
sion factor along sample lines-of-sight vs. A/comoving 
distance from the most massive galaxy in the compu- 
tational volume. There is significant transmission in 
the proximity zone and transmission gaps in the mean 
IGM away from the source. Gas infall results in some 
absorption red- ward of the line center. 



The observations of hig h-redshift QSO's ('Fan et al 



20011 : iBecker et al.l l200lh an d large-scale CMB polar- 



ization (|Spergel et al.l 120071 ) indicate that the inter- 



galactic medium has been completely ionized by red- 
shift z Q through an extended process. The most 
probable cause was the ionizing radiation of the First 
Stars and QSO's. Currently these are the two main 
direct observational constraints on this epoch. This 
scarcity of observational data is set to change dra- 
matically in the next few years, however. A num- 
ber of large observational projects are currently un- 



of hydrog;en (e.g . 



Madau et al. 


1997; 


Tozzi et al. 2000: 


Mellema et al. 


2006: iFurlanetto et al. 



2006), detection of small-scale CMB anisotropies due 
to the kinetic Sunyaev-Zerdovich (kSZ) e ffect (e.g. iHu 
boOOHSantos et ahbood : llliev et al.ll2007ch . and surveys 
of high-redshift Ly-g emitters and studies of the IGM 



absorption fe.g. | RhQads et al.ll2003l : IStanwav and et al 



20041 : iBunker et al 



The planning and success 
of these experiments relies critically upon understand- 
ing the large-scale geometry of reionization, i.e. the 
size- and spatial distribution of the ionized and neu- 
tral patches. This is best derived by large-scale sim- 
ulations, although a nu mber of semi- analytica l mod- 



els exist as well (e.g. IFurlanetto et al.l [2004). Re- 



cently we presented the first large-scale, high-resolution 
radiative t ransfe r simulations of cosmic reionization 
dlliev et al.ll2006l . [2007bh and applied those to derive a 



2006; 



range of reionization observables ("Mellema et cdj 
Iliev et al 2007c : Holder et al. 2007; Pore et al 
Iliev et al.ll2007a| |d|). Here we summarize recent results 



20071; 



on the characteristic scales and topology of reionization 
and implications of our simulations for the observability 
of high-redshift Ly-a sources. 



2 Characteristic scales of Reionization 

The characteristic scales of the ionization and density 
fiuctuations are directly imprinted in the fiuctuations of 
all reionization observables (Fig.[T]). It is thus critical to 
understand what these scales are and their dependences 
on the (still not yet well-known) reionization parame- 
ters. We investigated the dependence of the character- 
istic sizes of ionized and neutral patches on some basic 
properties of the sources of reionization and the inter- 
galactic medium, such as the halo mass-to-light ratio, 
susceptibility of haloes to positive and negative feed- 
back, and the gas clumping at small scales. We used 
two independent methods for identifying the size distri- 
bution of patches, the friends-of-friends (FOF) method 
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Fig. 5. — Ly-a sources at redshift z = 6.0: Luminosity 
function without (black) and with absorption included 
(red). For reference, the green, dotted line shows the 
result suppressing each source by 50%, which would be 
the case if e.g. all of the blue wing of the emission line, 
but none of the red wing is absorbed. 



( Iliev et al.l |2006[ ). where all topologically-connected 
ionized regions are consid ered as one re^o n, and the 
spherical average method (|Zahn et al.ll2007[ ). where an 
averaging over a spherical region and an ionization 
threshold are used to define the size. In the FOF 
method, throughout most of the evolution there is one 
very large connected region in which most of the vol- 
ume is contained. For the spherical average method, 
the bubble distribution typically peaks on Mpc scales. 
Suppression of ionizing sources within ionized regions 
and gas clumping both reduce the size and increase the 
number of H II regions, although the effect is modest, 
reducing the typical radius of H II regions by factors of 
at most a few (Fig. [2]). We also found that density and 
ionized fraction are correlated on large scales, regardless 
of the degree of clumping and suppression. The genus of 
the ionization field proved to be much more sensitive to 
suppression and clumping than the size distributions or 
the power spectra and could thus be used for more de- 
tailed characterization of the reionization parameters. 
We will present these results in detail in an upcoming 
paper (Alvarez et al., in prep.). 



3 High-redshift Ly-a sources 

Observations of the high-redshift Ly-a sources have 
provided us with a wealth of information about the 
state of the IGM and the nature of the galaxies at 
the end of reionization and hold even more promise for 



the future. In addition to probing the IGM and the 
source luminosity function (and thus, indirectly, the 
halo mass function and ionization source properties), 
they may also be used to constrain the reionization 
topology (jRhoadsl 120071 ). The first rare objects form 
in the highest peaks of the density field. The statistics 
of Gaussian fields predicts that such peaks are strongly 
clustered at high redshift. Hence each high-redshift, 
massive galaxy was surrounded by numerous smaller 
ionizing sources. The self-consistent simulations of such 
regions require following a large volume while also re- 
solving the low-mass halos driving reionization. 

In Figure [3] we illustrate several stages of the reion- 
ization history of a high density peak. The most mas- 
sive source in our volume (at z = 6) is shifted to the 
centre using the periodicity of the computational box. 
At redshift z — 12.9 the source is invisible due to the 
damping wing from the neutral gas outside the small 
H II region. By redshift z = 9 many more haloes 
have formed, most of them in large clustered groups. 
The H II region surrounding the central peak is among 
the largest, but the central source emission remains 
strongly affected by damping. Only by redshift z = 8 
the ionized region is large enough to render the source 
potentially visible. The reionization geometry becomes 
quite complex, most ionized bubbles are interconnected, 
but large neutral patches remain between them, which 
is also showing in the topological characteristics of the 
filed like the genus discussed above. Finally, by the 
nominal overlap z = 6.6 all ionized regions have merged 
into one topologically-connected region, although sub- 
stantial neutral patches remain interspersed throughout 
our volume, which is still largely optically-thick to both 
Ly-a and ionizing radiation. Only by z ~ 6 this volume 
becomes on average optically-thin to ionizing radiation. 

In Figure [H we show some sample absorption spec- 
tra for the same luminous source. The spectra exhibit 
extended high-transmission (10-60% transmission) re- 
gions in the highly-ionized proximity zone of the lu- 
minous source, within 5 Mpc/i~^ (~ 20 A). The cen- 
ter of the peak itself is optically-thick due to its high 
density. The infall around the central peak blue-shifts 
photons, resulting in some transmission behind the 
redshift-space position of the source. Away from the 
proximity region the absorption is largely saturated, 
but there are a number of transmission gaps with up 
to a few per cent transmission. In Figure [5] we show 
the Ly-a source luminosity function at the same red- 
shift. For the weaker sources roughly half of the intrin- 
sic luminosity is transmitted (the red wing of the line), 
while the most luminous sources suffer from additional 
absorption due to the gas infall that surrounds them. 
Future work will quantify the statistics of these features 
and its evolution. 
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